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In this paper we discuss two mechanisms by which high energy electrons resulting from dark 
matter annihilations in or near the Sun can arrive at the Earth. Specifically, electrons can escape 
the sun if DM annihilates into long-lived states, or if dark matter scatters inelastically, which would 
leave a halo of dark matter outside of the sun. Such a localized source of electrons may affect 
the spectra observed by experiments with narrower fields of view oriented towards the sun, such 
as ATIC, differently from those with larger fields of view such as Fermi. We suggest a simple test 
of these possibilities with existing Fermi data that is more sensitive than limits from final state 
radiation. If observed, such a signal will constitute an unequivocal signature of dark matter. 
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I. DARK MATTER IN THE SUN 

High-energy particles from dark matter (DM) capture 
and annihilation in the Sun offer a striking signature of 
dark matter [TJ [2] . The study of energetic neutrinos from 
the Sun [3H5] has received great attention in this context, 
as it is assumed that charged products would not escape 
the Sun's interior. Recent data and theoretical develop- 
ments call this assumption into question. In particular, 
the solar signatures of dark matter annihilation in the 
Sun can be greatly altered for dark matter that annihi- 
lates into a new force carrier [SHE], or for inelastically 
interacting dark matter (iDM) [S]. In this paper, we dis- 
cuss how cither scenario allows charged particles from 
DM annihilations in the Sun to reach the Earth, and the 
observational signatures of this effect. 

In the first case, illustrated in Fig. [ija), DM annihi- 
lates into long-lived particles, such as scalars associated 
with a new gauge sector. These long-lived particles can 
easily escape the Sun, and their subsequent decay in the 
solar system into electrons, muons, or charged pions can 
be detected. In the second case, DM captured through 
inelastic scattering may lack the minimum kinetic energy 
required to scatter again. If the elastic scattering cross- 
section is small, DM forms a loosely bound halo around 
the Sun and can annihilate outside the Sun as shown in 
Fig. [3b). 

In cither scenario, satellite observatories such as 
Fermi [TU] can detect the electronic annihilation prod- 
ucts as a cosmic ray electron excess strongly correlated 
with the Sun's direction. If observed, such an effect is an 
unequivocal signature of DM since no known astrophysi- 
cal phenomena can generate such a high-energy electron 
flux from the Sun. This type of signature may offer a 
unique probe of inelastically interacting dark matter, for 
which direct detection constraints are quite weak. 

Our estimates will show that a solar flux F ~ 
10~ 4 m _2 s _1 of particles above several hundred GeV 
should be detectable by experiments such as Fermi. 




FIG. 1: Two possible escape mechanisms for high energy 
charged particles from DM annihilations in the Sun. (a) DM 
may annihilate into long-lived states which first escape the 
Sun and only later decay, (b) DM may annihilate outside the 
Sun. 

Thus, only a small fraction of DM captured in the sun 
must annihilate through these channels to observe an ef- 
fect. Indeed, if for a given DM mass we take the largest 
cross-section allowed by direct detection limits on spin- 
independent elastic scattering (a$i ss 0.5 (3) x 10 -43 cm 2 
for m x w 0.1 (1) TeV) [TO [12], then DM is captured at a 
rate [IS], 

C «1.4x 10 21 s- 1 (^) 2/3 - (1) 

IDM models allow much larger cross-sections a n > 
10~ 40 cm 2 and hence considerably higher capture 
rates [HI [15]. For cross-sections in this range, the DM 
density accumulated over the age of the Sun is high 
enough that DM capture and annihilation rate (Ta) 
reach equilibrium so that Ta — \Cq. Assuming one 
observable product per annihilation actually leaves the 
Sun, the flux at the Earth is 

F ~ 5 x 10~ 3 urV 1 (elastic) (2) 
F ~ 50 m" 2 s" 1 (inelastic). (3) 

Both estimates are significantly larger than the sensitiv- 
ity limit, so that even very sub-dominant reactions of 
the form discussed here can be probed by careful analy- 
sis of electronic cosmic ray data. Interestingly, a signal as 
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FIG. 2: Geometry of DM annihilations into long-lived states, 
followed by 2-body decays into charged particles. Oi a b is the 
laboratory angle of the forward going decay product. 



large as ^ would constitute a significant fraction ~ 10% 
of the electronic cosmic rays above 300 GeV observed by 
PPB-BETS p3], ATIC [IZ|, and Fermi, and could lead 
to an observable feature in the overall signal rate. This 
is particularly important when considering balloon-based 
experiments situated at the south pole, such as ATIC and 
PPB-BETS, because spectral sculpting can be large for 
signals with a sharp directionality. 



II. ESCAPE MECHANISMS FOR CHARGED 
PARTICLES FROM SOLAR DM ANNIHILATION 

A. Long-Lived States 

We first consider the case as initially proposed in [5] 
where DM annihilates into the SM via some intermedi- 
ate states that then decay to Standard Model particles 
through a small mixing. This is the case, for example, if 
DM interacts through some new force, as has been em- 
ployed to explain astrophysical anomalies [BHEl 1X51 - 120] . 
These new bosons can decay back into the SM or cas- 
cade further in the dark sector. Some of the interme- 
diate states may be sufficiently long-lived that they es- 
cape the Sun before decaying into SM states [2TH23] (see 
[2HHS] for existing constraints on such long-lived states). 
As suggested in [37J , Lithium abundance discrepan- 
cies [28H30] may be resolved by late-decaying (r ~ 10 3 s) 
particles of weak-scale mass. This offers another moti- 
vation, as DM annihilations into such an intermediate 
state would allow the eventual decay products to escape 
the Sun. 

Fig. [2] illustrates the geometry of our setup and defines 
notation. To keep the discussion somewhat general, we 
calculate the profile of decay products assuming that DM 
annihilates in the center of the Sun into a particle with 
proper lifetime t, mass to/, and velocity /?. For simplic- 
ity, we assume that the intermediate particle decays via 
a 2-body reaction into an electron-positron pair. From 
Fig. [2j we see that the rate per unit area, per unit solid 
angle, per unit energy of electrons observed by a detector 
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where 9 cm and 9' are related via boosting from the CM 



frame to the solar frame, 
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distribution of the intermediate state given in its rest 
frame, and all quantities are evaluated in the geometry 
of Fig. [2] The energy delta function is simply enforcing 
that the solar frame energy is the energy of the boosted 
decay products. In particular, 

E(9, D, R) = imi ^{l + (3 cos 9 cm (9,R,D)). (5) 

Integrating over R and 4> can be done using the delta 
function and assuming azimuthal symmetry. We treat 
7, P and D as fixed. The delta function localizes R as 
a function of E (where E is the energy of the observed 
lepton), and 9. Also note that on the support of the delta 
function, 9' is a function of E alone, and r is a function of 
E and 9. Introducing the dimensionless energy variable, 
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we have that the minimum and maximum value for x is, 
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where the last equality follows simply from the law of 
sines: r=^|*. 

Putting all the pieces together, and using d £f g — = 
1/2 for isotropic 2-body decay, we have, 
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where x max 
L = jcPt. 

For P « 1, and 7CT < AU, the decay products travel 
in approximately straight lines from the Sun, and so the 
profile is obviously peaked strongly at the Solar center 
relative to the detector. For P <1 and jcPt > AU, the 
directionality relative to the Sun is significantly broad- 
ened — see Fig. [3] 
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dN/dfl : y=20 and 1.25, cr=.0.1. .and 5 AU 
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FIG. 3: Angular distribution from a single 2-body cascade, 
in degrees, for various decay lengths and long-lived particle 
velocities. We consider 7 = 1.25 (red) and 20 (blue), and 
jct = 0.1 (dashed) and 5 (solid) AU. If dark matter anni- 
hilates through a new force at the GeV-scale, even larger 7 
are expected. Both decay length > 1 AU and moderate 7 
are required to broaden the electron and FSR photon fluxes 
significantly. 



B. DM Annihilations Outside the Sun 

Another mechanism allowing electronic annihilation 
products to reach the Earth is for DM to annihilate 
near the surface of the Sun. In the standard scenario 
of elastic scattering, DM undergoes many further colli- 
sions with matter in the Sun after being captured. If the 
spin-independent elastic scattering cross-section is larger 
than about 10 -47 cm 2 the collision rate is sufficiently 
high that the captured DM will quickly thermalize with 
the rest of the matter in the Sun and concentrate in the 
inner core [5] (this estimate is based on the iron density 
and abundance in the sun). We therefore typically do 
not expect appreciable annihilation rate outside the Sun, 
as was clearly demonstrated in Ref. [31] . 

The situation is quite different in the iDM scenario. In 
particular, DM particles scatter only a few times before 
they lose enough energy to render further scattering kinc- 
matically forbidden. If the elastic scattering component 
is much smaller than 10~ 47 cm 2 , the DM will never reach 
thermal equilibrium with matter in the Sun resulting in 
only mild core concentration [33]. A non- negligible frac- 
tion of captured DM is then bound in elliptical orbits of 
order the size of the Sun and can then annihilate outside 
of the Sun. This is different from the usual case of elastic 
scattering, which was carefully investigated in Ref. [3T] 
in that the WIMPs can spend a very long (of the order of 
the lifetime of the Sun) time outside the Sun since they 
never lose enough energy to fall into the center and get 
trapped. 

The iDM capture rate in the Sun was calculated in 
[14l 115] . While the DM no longer thermalizes with mat- 
ter in the Sun, the annihilation rate does track the cap- 
ture rate since the WIMPs will continue to accumulate 
until equilibrium is reached. Using er„ = 10 -40 cm 2 , the 
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FIG. 4: A plot of / ou t against the inelasticity parameter 5 
for m x = TeV. The dependence on m x is weak since the 
dynamics is determined only by the reduced mass of the DM 
particle-nucleus system. For much higher values of the inelas- 
ticity, the capture rate rapidly drops as it becomes impossible 
for the DM particles to scatter anywhere in the Sun. 



annihilation rate outside the Sun is then related to the 
capture rate as, 
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where a conservative estimate [44] of / ou t is, 
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The resulting flux of charged particles at the Earth is, 
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To compute / ou t (shown in Fig. |4j) , we simulated the 
accumulation of DM in the Sun assuming only inelastic 
collisions as described in Ref. [13]. Notice that even at 
low inelastic thresholds with / out ~ 10 -4 , the annihila- 



tion rate outside the Sun (eq. 12 ) and resulting charged 
particle flux (eq. 14) can be quite large — much larger 



than our estimated Fermi's sensitivity. 



III. EXISTING CONSTRAINTS AND 
SIGNATURES 

A. FSR Constraints 

Any electronic production near the Sun or along our 
line of sight also yields 7-rays from final state radia- 
tion (FSR). Therefore, these processes are bounded by 
Fermi's measurement of the 7-ray spectrum from the 
Sun up to ~ lOGeV [33]. For illustration we consider 
DM annihilation XX ~^ 02 02 with two possible decays 
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for <\>2- one-step decay (f>2 —> e + e~; and two-step de- 
cay 4>2 — ► 4>i4>i followed by <pi — ¥ e + e~. The limits 
have been obtained in [MJ US] [33] using the simplified 
expressions for differential photon flux from |36j . Com- 
parable limits can be obtained from the Milagro search 
for high-energy gamma rays, which is sensitive to very 
high-energy (> 100 GeV) [37] ■ These constraints con- 
sequently severely restrict high-energy electronic emis- 
sions from the Sun. For example, an injection spectrum 
normalized to the positron flux observed by PAMELA 
would produce a 7-ray flux at 5-10 GeV two orders of 
magnitude larger than observed. It is therefore difficult 
to explain the low energy positronic excess as coming 
from annihilations of DM in the Sun without running 
into a serious conflict with the Fermi 7-ray observations. 
Likewise, attempting to explain the discrepancy between 
ATIC and Fermi data by normalizing the electron flux to 
the difference in flux at higher energies E ~ 500 GeV, the 
resulting 7-ray flux is still about an order of magnitude 
too large. FSR constraints generically restrict the Solar 
electronic flux to be at most O(10%) of the total flux ob- 
served by Fermi. Such a rate is still significantly above 
the expected sensitivity of a dedicated search for high- 
energy electrons from the Sun. Moreover, much larger 
electronic fluxes are possible if m < j >1 « 2m e , in which case 
FSR is phase-space suppressed, or if the injection profile 
is spatially broad, as is the case for /3 < 1, 7CT > AU (see 
Fig.§. 

B. Electronic Signatures 

Testing the scenarios described above is possible with 
existing Fermi electron data [TU]. One possibility is to 
form an asymmetry variable based on the differential flux 
coming from outside of the sphere defined by the Sun 
(night) or from the inside (day). One can then eliminate 
much of the systematics and be sensitive to any residual 
flux from the general direction of the Sun. 

A dn = d$ e /dE\ day - d<P e /dE\ night (15) 

Before discussing some more precise signatures, we briefly 
consider the trajectory of an electron emitted from the 
sun. The strong magnetic fields in the Sun's vicinity [38] 
will affect any electronic flux emanating from the Sun. 
Since the gyromagnetic radius is smaller than the curva- 
ture in the fields the electrons will follow the field lines, 
causing them to diffuse and deflect. Very close to the 
Sun's surface, an approximate dipole field of 1 gauss is 
present which will result in a diffusion of the isotropic 
flux. Away from the Sun's surface, the well known Parker 
spirals 39J will also alter the electrons' trajectories caus- 
ing them to arrive at the Earth from an apparent source 
shifted by up to 30° relative to the Sun's position. These 
two effects mean that we must consider a few possibilities 
when searching for electrons from the sun. In the pres- 
ence of a positive signal, these effects will clearly need to 
be better understood, as well as the effects of the Earth's 



magnetic field on the acceptance of experiments situated 
at the poles, such as ATIC and PPB-BETS, in order to 
interpret any result. 

However, there are simple and straightforward ap- 
proaches that can increase sensitivity without a detailed 
modeling of the signal. In particular, we can con- 
sider the more precise signature of the differential flux, 
d$sun/dE(A8) in a circle of size AO centered on the Sun. 
By comparing to a region of the sky of the same size but 
an angular distance away (a "fake Sun") , an estimate 
of the expected background flux can be obtained. (Opti- 
mally the fake sun should be in the hemisphere opposite 
that of the real sun.) If <J>f a ke is the flux through the "fake 
Sun" , then a useful search variable is the asymmetry, 

Ag = d$su tt (A0)/dE - d$ {3ke {A6)/dE (16) 

studied as a function of energy. We expect that the 
background rates will be well approximated by the "fake 
Sun" for regions of size AO > 10 degrees at energies above 
~ 100 GeV. For smaller circles and lower energies, So- 
lar magnetic fields and Earth magnetic field asymmetries 
may alter the expectation beyond statistical uncertain- 
ties. In Fig. [5] we plot the differential flux for several 
benchmark signals, namely one- and two-step cascade 
decays for ITeV and 2TeV DM mass. The fluxes are 
normalized so that the FSR differential photon flux at 
10 GeV is d^-y/dE = 3 x 10- 6 /m 2 /s/GeV, hence satu- 
rating the last bin of the Fermi solar result [34 . The 
fluxes for the different benchmarks are given in Table [T] 
Fig. [5] also includes background estimates with statistical 
uncertainty for a few choices of AO. Given the figures for 
the different models, Solar electronic fluxes as small as 
F ~ 10~ 4 m^s" 1 at energies E> 100 GeV should be 
detectable with existing Fermi data. 

TABLE I: Fluxes of the annihilation products from the Sun 
used in Fig. [5] The fluxes are normalized so that the 
FSR differential photon flux at 10 GeV is d$ y /dE = 3 x 
10 _6 /m 2 /s/GeV saturating the last bin of the Fermi solar re- 
sult [33]. To obtain the actual electronic flux these figures 
should be multiplied by 2 (4) for the one- (two-) step cas- 
cades, assuming the annihilation products decay entirely into 
e + e~ pairs. 



Model 


Flux (m- 2 s _1 ) 


one-step, 1 TeV 
one-step 2 TeV 
two-step 1 TeV 
two-step 2 TeV 


2.7 x 1(T 4 

2.5 x 1(T 4 

1.6 x 1(T 4 
1.4 x 1(T 4 



Such an approach helps with the defocusing effects of 
the Solar magnetic fields on the electron fluxes, but what 
about the possibilities of a coherent deflection from the 
Parker spirals? A simple approach for this would be to 
look for "hotspots" in the solar frame. That is, one can 
consider the flux in regions of size 6 which tile the 

sky. In observing this distribution, there will certainly 
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FIG. 5: Electronic energy spectrum for a one-step (dashed) 
and two-step (solid) decays for a DM mass of 1 TeV (red) 
and 2 TeV (blue). These correspond to the process \X ~^ 
4>2<l ) 2 followed by (j>2 — > e + e~ (one-step) or <p2 ^ 4>i4>i an d 
(pi — > e + e~ (two-step). The curves have been normalized to 
avoid FSR limits as discussed in the text. Also plotted are 
estimates of Fermi's statistical sensitivity for half sky (upper) 
and 30° searches, given the statistics of [10] 



regions around the sun, one should be careful in placing 
limits from looking at regions smaller than about 30°. 
Should the signal arrive from a different angle, highly 
focused regions could give anomalously strong limits on 
electronic production. 



IV. CONCLUSIONS 

In this paper we proposed two novel mechanisms by 
which high energy electrons originating from DM anni- 
hilations in the Sun can arrive at the Earth's vicinity. 
Such high energy electron flux should be observable al- 
ready with existing Fermi data and may help resolve the 
outstanding discrepancy with the ATIC results. If such 
a flux is observed, it will constitute an unambiguous evi- 
dence for DM annihilations in the Sun. It will also force 
a revision in our ideas about the way DM interacts with 
matter and/or the nature of its annihilation products. 



be variation from place to place, but the appearance of 
outliers (i.e., anomalously large d ^ d ^ in one direction 
compared to the others) could give evidence of solar pro- 
duction, even if this hotspot were not centered on the 
sun. The expectation is that while some deflection is 
possible, extremely large deflections ( > 30°) are unlikely. 
Although this technique is less sensitive, it overcomes the 
complications associated with a coherent deflection of the 
electronic flux. As a final corollary to this point, we note 
that while discoveries can be made by looking at smaller 



Acknowledgments 

We thank J. Wefel for useful discussions of the ATIC 
experiment, and S. Murgia and N. Giglietto for useful 
discussions of the Fermi solar gamma-ray measurements. 
PS is supported by the US DOE under contract number 
DE-AC02-76SF00515. NW is supported by NSF CA- 
REER grant PHY-0449818 and DOE OJI grant # DE- 
FG02-06ER41417. IY is supported by the James Arthur 
Fellowship. 



W. H. Press and D. N. Spergel, Astrophys. J., 296, 679 [13 
(1985). 

J. Silk, K. A. Olive, and M. Srednicki, Phys. Rev. Lett., 
55, 257 (1985). [14 
L. M. Krauss, M. Srednicki, and F. Wilczek, Phys. Rev., 
D33, 2079 (1986). [15 
T. K. Gaisser, G. Steigman, and S. Tilav, Phys. Rev., 
D34, 2206 (1986). [16 
K. Griest and D. Seckel, Nucl. Phys., B283, 681 (1987). 
D. P. F inkbeiner and N. Weiner , Phys. Rev., D76, 083519 [17 

(2007) , [arXrv:astro-ph/0702587| . ' [18 
M. Pospelov, A. Ritz, and M. B. Voloshin, Phys. Lett., 
B662, 53 (2008), |arXiv:0711.4866 [hep-ph]| . [19 
N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, 

and N. Weiner, Phys. Rev., D79, 015014 (2009), [20 
|arXiv:0810 0713 [hep-ph]] . 

D. Tucker-Sm ith and N. Weiner, P hys. Rev., D64, [21 
043502 (2001), [arXiv:hep-ph/0101138| . 
A. A. Abdo et al. ( The Fermi LATJ7 Phys. Rev. L ett., 
102, 181101 (2009), | arXiv:0905.0025 [astro-ph.HE]] . [22 
J. Angle et al. (XENON) , Phys. Re v. Lett., 100, 021303 

(2008) , [arXiv:0706. 0039 [astro-ph]| . [23 

Z. Ahmed et al. (CDMS), (2008), |arXiv:0802. 3530 [astro-| 

"phU . [24 



This formula gives an excellent approximation to the cap- 
ture obtained from the more precise computation found 
in gDHH]. 

S. Nussinov, L. T. Wa ng, and I. Yavin, (2009), 
|arXiv:0905.1333 [hep-ph]] . 

A. Menon, R. Morris, A. P ierce, and N. Weiner, (2009) , 
|arXiv:0905.1847 [hep-ph]] . 

S. Torii e t al. (PPB-BETS), (2008), |arXiv:0809.0760] 
[astro-ph] . 

J. Chang et al, Nature, 456, 362 (2008). 

I. Cholis, L. Goodenough , and N. Weiner, (2008), 

|arXiv:0802.2922 [astro^phj] . 

M. Pospelov and A. Ritz Phys. Lett., B671, 391 (2009), 
|arXiv:0810.1502 [hep-ph]] . 

I. Cholis, D. P. Finkbeiner, L. Goodenough, and 
N. Weiner, (2008), [arXiv:0810. 5344 [astro-ph]| . 
M. Baumgart, C. Cheung, J. T. Ruderman, L.-T. Wang, 
and I. Y avin, JHEP, 04, 014 (2009), |arXiv:0901. 0283 
[hep-ph] | . 

B. Batell, M. P ospelov 
|arXiv:0903.0363 [hep-ph]] . 
R. Essig, P. Schuster, and N 



and A. Ritz, (2009), 
Toro, (2009), 



arXiv:0903.3941 [hep-ph] . 
P. Schuster, N. Toro, and I. Yavin, Phys. Rev. 



D81, 



G 



016002 (2010), |arXiv:0910. 1602 [hep-ph]| . 
[25] J. D. Bjor ken, R. Essig, P. Schuster, and N. Toro, 

(2009) , arXiv:0906.0580 [hep-ph] . 

[26] K. Jedamzik, Phys. R ev. Lett., 84, 3248 (2000), 

|arXiv:astro-ph/9909445| . 
[27] K. Jedamzik, Phys. R ev., D77, 063524 (2008), 

|arXiv:0707.2070 [astro-ph]| . 
[28] E. Vangioni-Flam et a/., (1998), arXiv:astro-ph/9811327 

[29] M. Asplund, D. L. Lambert, P. E. Nissen, F. Primas, and 
V. V. Smith Astrophys. J., 644, 229 (2006), |arXiv:astro- 
ph/0510636] . 

[30] R. H. Cyburt, B P. Fields, and K. A. Ol ive, JCAP, 

0811, 012 (2008), |arXiv:0808.2818 [astro-ph]| . 
[31] S. Sive rtsson and J. Edsjo, Phys. Re v., D81, 063502 

(2010) , |arXiv:0910. 0017 [astro-ph.HE]] . 

[32] Such small elastic cross-sections are expected in cases 
where the inelasticity is achieved by lifting the degener- 
acy between the two CP eigenstates of a scalar or the two 
Majorana components of a Dirac fermion. 

[33] The quantity / out was calculated based on the calcula- 
tions presented in Ref. [T2]- An analytical solution for 
the trajectories inside and outside the Sun was used in a 
Monte-Carlo simulation of the capture and entrapment 
process. The simulation results, which include the final 
trajectory of a large number of particles, were used to 
calculate the asymptotic density for a give WIMP mass, 
m x and inelasticity 5. 



[34] N. Giglietto, TeV Particle Astrophysics Conf. (2009). 
[35] B. Batell, M. Pospelov, A . Ritz, and Y. Shang, (2009), 

|arXiv:0910.1567 [hep-ph]] . 
[36] J. Mardon, Y. Nomura, D. Stolarski, and J. Thaler, 

JCAP, 0905, 016 (2009), [arXiv:0901. 2926 [hep-ph]| . 
[37] R. W. Atkins et al, Phys. Rev., D70, 083516 (2004). 
[38] H. W. Babcock and H. D. Babcock, Astrophys. J. , 121, 

349 (1955). 

[39] E. N. Parker, Astrophys. J. , 128, 664 (1958). 
[40] A. Gould, Astrophys. J., 321, 571 (1987). 
[41] P. Gon dolo et al, JCAP, 0407, 008 (2004), |arXiv:astro- 
ph/0406204| . 

[42] This formula gives an excellent approximation to the cap- 
ture obtained from the more precise computation found 
in HUE]. 

[43] Such small elastic cross-sections are expected in cases 
where the inelasticity is achieved by lifting the degener- 
acy between the two CP eigenstates of a scalar or the two 
Majorana components of a Dirac fermion. 

[44] The quantity / out was calculated based on the calcula- 
tions presented in Ref. [14]. An analytical solution for 
the trajectories inside and outside the Sun was used in a 
Monte-Carlo simulation of the capture and entrapment 
process. The simulation results, which include the final 
trajectory of a large number of particles, were used to 
calculate the asymptotic density for a give WIMP mass, 
m x and inelasticity 5. 



